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Inspired by a recent experiment of localization-delocalization transition (LDT) of indirect excitons 
in lateral electrostatic lattices [M. Remeika et al, Phys. Rev. Lett. 102, 186803 (2009)], we 
investigate the interplay between periodic potential and nonlinear interactions of indirect excitons 
in coupled quantum wells. It is shown that the model involving both attractive two-body and 
repulsive three-body interactions can lead to a natural account for the LDT of excitons across the 
lattice when reducing lattice amplitude or increasing particle density. In addition, the observations 
that the smooth component of the photoluminescent energy increases with increasing exciton density 
and exciton interaction energy is close to the lattice amplitude at the transition are also qualitatively 
explained. Our model provides an alternative way for understanding the underlying physics of the 
exciton dynamics in lattice potential wells. 

PACS numbers: 73.63.Hs, 78.67.De, 71.35.Lk, 73.20.Mf 



I. INTRODUCTION 

Indirect excitons (spatially separated electron-hole 
pairs) in coupled quantum wells (CQWs) can have a long 
lifetime and a high cooling rate. With these two merits, 
Butov et al. have successfully cooled the trapped ex- 
citons to the order of IK and have observed surprisedly 
that excitons form two, inner and external, rings to which 
periodic bright spots appear in the external ring.^'^ Al- 
though there is no clear evidence that these excitons are 
condensed into the BEG state, it is fascinating enough to 
see some puzzling particle number distributions in var- 
ious confined potential wells<^i^ More recently, periodic 
potential (lattices) due to gate voltages was created for 
indirect excitons of built-in electric dipole moment. It 
enables observation of localization-delocalization transi- 
tion (LDT) for transport across the lattice with reduc- 
ing lattice amplitude or increasing particle density.^ This 
gives an opportunity to understand further the complex 
dynamical behaviors of indirect excitons. 

In the literature, a charge separated transportation 
mechanism was proposed,— which gave a satisfactory ex- 
planation to the formation of exciton rings and the dark 
region between the inner and the external rings. How- 
ever the origin of the periodic bright spots in the external 
ring can not be fully understood within the above frame- 
work. Alternatively a self-trapped interaction model in- 
volving attractive two-body and repulsive three-body in- 
teractions in the system was proposed ji^ It was shown 
that the interplay between the two-body attraction and 
the three-body repulsion can give a good account to the 
periodic bright spots in the external ring. In addition, 
the self-trapped interaction model also explained well 
the abnormal exciton distribution in an impurity poten- 
tial, where the photoluminescence (PL) pattern becomes 



much more compact than a Gaussian with a central in- 
tensity dip, exhibiting an annular shape with a darker 
central region^S Moreover, the model also captured some 
other experimental details, for instances, the dip can turn 
into a tip at the center of the annular cloud when the 
sample is excited by higher power lasers. 

The success of the self-trapped interaction model in 
understanding various phenomena has motivated us to 
investigate the LDT phenomenon reported in Ref. [5| us- 
ing the same model. It will be shown that the complex 
LDT as well as other dynamic behaviors can also be qual- 
itatively explained by the self-trapped interaction model. 

This paper is organized as follows. In Sec. |lll we in- 
troduce the self-trapped interaction model described by 
a phenomenological nonlinear Schrodingcr equation. An 
important factor regarding how the exciton distribution 
depends on temperatures and energies is detailed. In 
Sec. mil numerical calculations and detailed discussions 
are given for studying the LDT of transports across the 
lattice. It will be shown that the model introduced in 
Sec. In] can lead to a reasonably good explanation for 
LDT. Sec. lIVl is a brief summary. 



II. MODEL 

The idea of the self-trapping model came from the den- 
sity dependence of various observed PL spectroscopy!^"— 
As a matter of fact, intensity of PL spectroscopy is pro- 
portional to the exciton number density and therefore 
complex experimental PL data will reveal complex inter- 
action between excitons. It is clear that the interaction 
between excitons is neither purely attractive, nor purely 
repulsive. If it is purely repulsive, it could drive the sys- 
tem towards homogeneous distribution. On the contrary. 



if it's purely attractive, the system will become unstable 
and eventually collapse when the exciton density is larger 
than some critical value. Thus it is reasonable to spec- 
ulate that at low densities the interaction is dominated 
by an attraction, while at high densities a repulsion must 
exist to prevent the system against collapses/ The ob- 
servation that exciton cloud first contracts then expands 
by increasing the excitation power gives a strong support 
to the above scenario.— 

Excitons behave as electric dipoles, so strong Coulomb 
repulsion will govern the indirect excitons in CQWs at 
high densities to which dipoles are aligned in paralleli^ 
However at low densities such that two dipoles change 
from aligning in parallel to inclined, the attraction be- 
tween the electron of one exciton and the hole of the other 
exciton will dominate insteadJ^ In addition to Coulomb 
interactions, exchange effect may be another important 
factor for interactions. When two excitons approach to 
each other, exchange interaction between two electrons 
(or two holes) will become more important. This may be 
another source of the attraction between excitons. 

In the dilute limit, it is assumed that indirect exci- 
tons in CQWs are governed by two-body attractions and 
three-body repulsions. Equivalently one can consider the 
following scaled nonlinear Schrodinger equatioi*^'^ 

-'^V^tpj{r) + [V,^{r)-gin + g2n^]tPj{r) ^ EjtPjir), (1) 

where tpj (r) and Ej are the j-th energy eigenfunction and 
eigenvalue and T4x denotes the external potential. The 
eigenfunction is normalized under J^ \ipj{r)\'^dQ, = 1. In 
Eq. ([1]) , all lengths and energies are scaled in units of ctpl 
and ?i^/m*cr|,L respectively, where m* = TOe/2 is effec- 
tive mass of the exciton and crpL is the root-mcan-square 
radius of the exciton cloud observed via PL.^ Consider- 
ing the special property of the system, local probability 
density of excitons [n = n{r)] will be taken to be 
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n(r)=^,7(i?,)|V,(r)| 



(2) 



where A^ denotes the total number of bound states in- 
volved and rjlEj) is the probability (or distribution) func- 
tion for energy-level Ej that satisfies the normalization 

condition J2i=i vi^j) = 1. At the mean- field level, two- 
and three-body coupling constants, gi and g2^ are defined 
positively, which will be proportional to N and iV^ with 
N being the number of excitons. Note that exciton num- 
ber density can be controlled by and proportional to the 
excitation laser power P, thus gi is proportional to P in 
the current model. 

Determination of the distribution function rj{Ej) in ^ 
requires considering both complex energy relaxation and 
recombination processes. At low exciton densities (n <C 
1/a^, aB being Bohr radius), relaxation due to exciton- 
exciton and exciton-carrier scattering can be neglected^i 
and consequently relaxation time is dominated by the 
scattering of excitons with acoustic phonons. At low bath 




FIG. 1. Calculated spatial probability density distribution 
of highly degenerate indirect excitons, n{x,y), for low den- 
sities {gi = 10) [frame (a) & (d)], intermediate densities 
{gi = 20) [frame (b) & (e)], and high densities {gi = 30) 
[frame (c) & (f)]. Frame (a)-(c) are with a lattice poten- 
tial Vex ~ uo cos(27ra;/A) {vo — 50 and A — 0.05), while (d)- 
(f) are without a lattice potential {vq = 0). In all frames, 
g2 = 0.001(71, a is fixed at 0.1, and x and y axes are in units 
of CTpL (see text) . The results are intended to be compared to 
those in Fig. 1 of Ref. [g. In frame (c) , FWHM (calculated in 
Fig. 3) and yo (used in Fig. 2) are denoted. 



temperatures, (Tf, < 1 K), this kind of relaxation rate 
decreases dramatically due to the low phonon densityJ^ 
For the recombination process, on the other hand, be- 
cause excitons in the lowest self-trapped level are quan- 
tum degenerate, they are dominated by the stimulated 
scattering when the occupation number is more than a 
critical value. Strong enhancement of the exciton scat- 
tering rate has been observed in the resonantly excited 
time-resolved PL experiment ji^ Therefore, even though 
the phonon scattering rate is larger than the radiative re- 
combination rate, it is in reality that the system does not 
reach thermal equilibrium. Following Ref. Q, in accord 
with the discrete PL spectroscopy reported in Ref. jl4| . 
ri{Ej) will have the following form 



v{Ej) ^Cexp[a{Ej ~ n)], 



(3) 



where C is the normalization factor, /x is the chemical 
potential, and a can be viewed as an "effective tempera- 
ture" . 



III. RESULTS AND DISCUSSIONS 



A. Spatial exciton distributions with or without 
lattices 



Key features of indirect excitons with or without lat- 
tices, as reported in Ref. p\, are the following. LDT 
for transport across the lattice was observed with reduc- 
ing lattice amplitude or increasing exciton density. The 
inner ring, which was observed in the PL patterns of 
indirect excitons without lattices, persists in the lattice 



case. For relatively low excitation power {P — 3.7 /iW), 
outer ring effect does not emerge because of low particle 
number density and hence low drift and diffusion of elec- 
trons and holes. However, for relatively higher excitation 
power {P = 12 ^W), in addition to the inner ring, outer 
ring seems to exist in the case without the lattices (see 
Fig. l-(iIof Ref. 0). With the lattices (see Fig. l.(f) 
of Ref. j5j ) , outer ring transforms to spatially separated 
clouds, while inner ring persists. 

Based on the above observations, in the current 
localized-delocalized system it seems that there are two 
kinds of excitons, whose spatial distributions could be 
governed by different mechanisms. The ones associated 
with the inner ring effect are ascribed to the hot excitons 
which have been studied in great details in Refs^^'^i^'i^. 
The other ones are ascribed to the cool excitons whose 
formation and pattern are considered to be due to quan- 
tum degeneracy.-'-. 

When excitation power P is low, densities of both hot 
and cool excitons are low. They are all localized within 
the range of the laser spot (i.e., localization). When P 
is increased, densities of both kinds of excitons will in- 
crease for which inner ring forms due to transportation 
and cooling of the hot excitons. While cool excitons will 
form patterns due to quantum degeneracy arising from 
complex interactions. Delocalization effect increases with 
increasing the laser power P, to which radius of the in- 
ner ring as well as the distance between the two sepa- 
rated degenerate exciton clouds increase. As a matter of 
fact, full width at half maximum (FWHM) of the exciton 
cloud will change with the particle number density or the 
lattice potential amplitude, as shown in Fig. 3 of Ref. [3- 

Here we shall focus on the LDT associated with highly 
degenerate low-energy (cool) excitons only. Spatial dis- 
tributions of low-energy exciton density are calculated 
self-consistently based on Eqs. (dJ-©. The results are 
shown in Figs. [T](a)-[Hc) for the case of a periodic lat- 
tice potential, Vc-x — vo cos(27ra::/A) with vq — 120 and 
A = 0.05, and in Figs. [lld)-[ljf) for the case without a 
lattice potential {vq = 0). In the case with lattice poten- 
tial, a weak rectangular potential well is applied to ensure 
the strip pattern along the lattices (i.e., with the Dirich- 
let boundary condition). In our model, highly degener- 
ate excitons are self-trapped due to competition between 
complex interactions and the nonequilibrium energy dis- 
tribution. Apart from the periodic lattice potential, any 
additional weak external potential will only have minor 
effect on the results. In our calculations, the wave func- 
tion is chosen to be a Gaussian initially. All self-trapped 
energy eigenstates ipj with Ej < are solved and will 
contribute to the spatial distribution. 

As mentioned before, at low excitation powers P low- 
energy exciton density profile essentially coincides with 
the excitation laser spot. In our model, this arises due to 
the dominant attraction between excitons in the dilute 
limit and at the same time only ground state contributes 
to the self-trapping. The ground-state wave function is 
s-wave with a maximum at the center [see Figs.[T{a) and 
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FIG. 2. (Color online) (a) Exciton probability density dis- 
tribution n{x,yo) with j/o denoted in Fig. [ijc). (b) Lattice 
potential energy density Ep{x,yo) = Vcx{x)n{x,yo) plotted 
vs. X. (c) Spatially dependent interaction energy density 
Ei{x,yo) and kinetic energy density Ek{x,yo) plotted vs. x. 
Parameters used are the same as those in Fig. [2c), except 
a larger gi = 50 and vq = 100 to account for the experi- 
ment. Vertical line shows that maxima in density n{x,yo) 
correspond to minima in periodic potential energy density 
Ep{x,yQ). In all frames, x axis is in units of ctpl. The results 
are intended to be compared to those in Fig. 2 of Ref. [g. 



[TJd)]. Excitons are localized and do not travel beyond 
the excitation spot. 

When exciton number is increased with higher excita- 
tion power P, excitons can delocalize and spread beyond 
the excitation spot [see Figs. [T](b), [T](c), [TJe), and (Hf)]. 
In these cases, effect of repulsion becomes more impor- 
tant which results more energy eigenstates involved in 
the self-trapping. The first excited states are nearly two- 
fold degenerate p-wave with a node at the center. Su- 
perposition of the ground-state s-wave and the two first- 
excited-state p-wave wavefunctions leads to an annular 
distribution with a hole at the center in the case without 
lattices [see Fig. [Ijf)]. With the lattices, two- fold degen- 
eracy is lifted and consequently ring-shaped delocalized 
excitons shift to two separate clouds at high excitation 
powers [see Fig. [T^c)]. 



B. Spatially dependent PL energy 

A PL image of indirect excitons in the delocalized 
regime, plotted as energy vs. space, is presented in 
Fig. 2(a) of Ref. [a. As a matter of fact, both integrated 
PL intensity I(x) and the average PL energy huj{x) show 
a small modulation at the lattice period superimposed 
on a smoothly varying profile. In particular, the average 
PL energy shows a dome carve shape. 

In fact in CQWs, spatially dependent energy of indi- 



rect excitons, which is directly related to the PL energy 
fiw, involves four major contributions. The first part is 
the intrinsic energy that includes contributions due to 
the band gap, the Coulomb interaction between electrons 
and holes, and the electric potential used to form indirect 
excitons. Intrinsic energy of exciton mainly contributes 
to the smooth part of the PL energy huj since it does 
not change with positions when experimental condition 
remains unchanged. The second part is the kinetic en- 
ergy Ek- Near the laser spot, excitons are hot with a 
large kinetic energy. It is believed that this is the main 
cause why a high energy distribution is located in the 
center regime. The third part is the interaction energy 
Ej. It also contributes to the smooth part of the PL 
energy huj. Within our model, when exciton interaction 
is in the repulsion-dominant regime, Ei increases with 
the exciton density n and thus smooth part of Huj will in- 
crease with increasing the excitation power P. Moreover, 
reducing the cloud size corresponds to increasing the exci- 
ton density when excitation power P remains unchanged. 
Therefore smooth part of hw will increase with reducing 
the cloud size. The last part is the lattice potential en- 
ergy Ep. Within the Thomas- Fermi approximation, it is 
easy to verify that Ep dominates over Ej and Ek in such 
a dilute system (see later). This is why maxima in exci- 
ton number density (and equivalently the PL intensity) 
correspond to minima in the periodic potential energy 
density Ep (see Fig. 2). 

Within the same model, average energy density and 
number distribution of highly degenerate excitons are cal- 
culated and shown in Fig. [2] These are intended to be 
compared to those reported in Fig. 2 of Ref. ^. Apart 
from the smooth part mainly due to the intrinsic energy, 
spatially dependent kinetic energy density of degenerate 
excitons can be calculated via* 
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while the interaction (mean-field) energy density is given 
by Ei{r) = —gin{r)^ + g2n{r)^ and the lattice potential 
energy density is given by Ep{r) = Vcx{j^)n{r) , which 
are obtained by multiplying ri{Ej)ip*{r) to Eq. ([T]) and 
summing over j. 

Fig. HJa) shows the calculated spatial probability den- 
sity distribution of degenerate excitons, n(a;, j/o)- Here 
yo corresponds to center of the upper cloud in y di- 
rection [see Fig. [D^c)]. The spatial lattice potential en- 
ergy density Ep{x,yo), kinetic energy density Ek{x,yo), 
and interaction energy density Ei{x, yo) are presented in 
Figs. [2{b) and [Ifc) respectively. In view of Fig. [2l Ek 
and Ei are about one order of magnitude smaller than 
Ep. This means that Thomas- Fermi approximation is 
valid in such a dilute exciton system. Consequently ex- 
citon number distribution is mainly determined by the 
periodic potential energy Ep. The higher the Ep is, the 
larger the local PL energy is. On the contrary, the lower 
the Ep is, the higher the local exciton number density is. 
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FIG. 3. (Color online) Calculated FWHM of delocalized exci- 
ton clouds vs. (a) logjQ(gi) ((j\ oc exciton number density n) 
for lattice amplitude vq = 50 and 100 and (b) wq for different n 
with g\ = 10 and 20. (c) Calculated number density n at the 
LDT point as a function of lattice amplitude dq. (d) Calcu- 
lated interaction energy density Ei{0, yo) at the LDT point as 
a function of lattice amplitude vq. {0,yo) corresponds to the 
outer cloud center. The results are intended to be compared 
to Fig. 3 in Ref. ^ 



As PL intensity is proportional to the local exciton num- 
ber density, it results that minima in energy correspond 
to the maxima in intensity. Equivalently the largest ex- 
citon number density (and hence the PL intensity) are 
located in the bottom of the periodic potential energy. 



Interplay between periodic potential and 
nonlinear interaction 



In view of the results shown in Figs. 3(a) and 3(b) of 
Ref. ^], it is surprising to see that FWHM of the exci- 
ton cloud (across the lattice) decreases with increasing 
the lattice amplitude. On the contrary, the excitation 
power needed for the transition from localized to delo- 
calized regimes increases with the increase of the lattice 
amplitudes [see Figs. 3(c) of Ref. i5i|. This is likely due 
to the reduction of the particle tunneling rate across the 
lattice. Thus a large laser power is needed to obtain ex- 
citon delocalized distribution when increasing the lattice 
amplitude. 

Fig. [3] studies the effects of number density and lattice 
amplitude on LDT of degenerate excitons. Here, as the 
dynamics of outer cloud is focused on, we define FWHM 
as the distance (the extension) between centers of the two 
clouds along the y-direction [see Fig. [TJc)]. In Figs.IDJa) 
and [S^b) , we show FWHM as a function of the particle 
number density (via changing gi) and the lattice poten- 
tial amplitude respectively. For weak excitation power 
P and low exciton number correspondingly, within our 
model both the attractive and repulsive interactions be- 
tween excitons are weak. All particles are sitting in the 



ground state and self-trapped in the shallow mean-field 
potential well. In this case, FWHM of exciton cloud is 
equal to zero and remains unchanged in the lower ex- 
citation power regime [see Fig. [Sja) and Fig. ^h) with 
5'i=10]. When the excitation power P or the particle 
density n is increased to above some critical point, LDT 
occurs. Within our model, repulsive interaction domi- 
nates and consequently excitons will be sitting in both 
ground state and two-fold degenerate first-excited states. 
Due to the one-dimensional lattice potential, two-fold de- 
generacy of excited states is lifted and consequently two 
separated delocalized exciton clouds form. 

In Fig. [SJc), we show the critical number density n 
for LDT as a function of lattice amplitude vq. When 
vq is increased, before the LDT transition excitons are 
strongly localized in a single well center, while particle 
tunneling rate across the well is largely decreased. As a 
consequence, critical number density for LDT increases. 
Within our model, before the LDT transition a larger 
lattice amplitude vq corresponds effectively (or relatively) 
to a smaller exciton density. This also means that the 
effect of the repulsion is decreased while the effect of the 
attraction is increased. This explains why a blue-shift of 
the LDT is observed in the experiment. It gives another 
evidence to support our model. 

In the delocalized regime, when excitation power P is 
low, only intrinsic energy and potential energy contribute 
to the exciton PL energy because interaction energy Ej 
is immaterial. When P is increases, however, Ej will 
become more important and contribute much more to 
the PL energy. Therefore, HAuj, defined as the difference 
between fku and that associated with the lowest P for 
LDT, is dominated by Ej. Within our model, in the de- 



localized repulsion-dominant regime, interaction energy 
density Ej will increase nionotonically with increasing 
the lattice amplitude vq. Fig. ^d) shows the calculated 
interaction energy density Ej{0,yo) at the LDT point as 
a function of lattice amplitude vq. Here {x,y) = (0,?/o) 
corresponds to the center of the outer cloud. In agree- 
ment with the experiment, it is seen that Ej at the LDT 
point is close to the value of the lattice amplitude vq. 
This agreement gives a strong support that the scaled 
parameters used in the current context are reasonable. 



IV. SUMMARY 

In summary, we have investigated the interplay be- 
tween periodic potential and nonlinear interaction of in- 
direct excitons in coupled quantum wells, with or without 
a lattice potential. The model which takes into account 
the competition between a two-body attraction and a 
three-body repulsion along with a reasonable nonequi- 
librium energy distribution gives an alternative quali- 
tatively good account on the localization-delocalization 
transitions (LDT) of excitons across the lattice when in- 
creasing the particle density or reducing the lattice am- 
plitude. 
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